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We have studied the fundamental process of hydrogen binding at interstitial, vacancy and grain 
boundary (GB) in palladium crystals using Density-Functional Theory. It showed that hydrogen 
prefers to occupy the octahedral interstitial site in Pd matrix, however a stable H-vacancy complex 
with most H occupations would contain up to eight hydrogen atoms surrounding the vacancy at 
tetrahedral sites. Furthermore, H presence assists the pairing or formation of nearby vacancies, which 
in agreement with previous suggestions by both experiment and theory investigation. Also, this 
observation could imply about a hydrogen embrittlement (HE) mechanism through the connections 
of microvoid and cracks. The segregation of hydrogen at grain boundary, nevertheless, has shown a 
different effect. High H accumulation results in grain boundary extension, which is related the HE 
mechanism of grain decohesion observed by experiments. 


I. INTRODUCTION 

The absorption of hydrogen in metals, in particular 
Pd-based materials, has attracted great interest due to 
its unique behavior, including its high solubility and ac- 
tivity 1-4 . However, it is also known that the interaction 
of impurity and the defects in host materials would dras¬ 
tically affect their properties and performance. Specif¬ 
ically, both experimental and theoretical investigations 
have evidenced that defects such as vacancy, dislocations 
and grain boundary in metals can attract a high local¬ 
ization of impurities, 5-9 and subsequently are responsible 
for many prominent phenomena. 

These effects, nevertheless, require a threshold concen¬ 
tration of hydrogen at intrinsic binding sites. Therefore, 
the learning of the H saturation at these locations could 
be critical for the understanding of their behavior. In 
general, hydrogen diffusion in crystals is active by first 
occupying interstitial positions until they are trapped in 
open volumes (such as vacancy or grain boundary). In 
fact, it was suggested that the coalescence of vacancies 
in the presence of H will provide source for microcrack 
initiation and subsequent embrittlement 10 . Also, exper¬ 
iments performed by Fukai et al. evidenced that high 
H absorption can introduce the formation of superabun¬ 
dant vacancies in several transition metals 6 ' 7 . Grain 
boundaries, on the other hand, can serve as ideal trans¬ 
port channels for impurity elements, at the same time 
could also be prefered nucleation site of dislocation and 
voids that results in stress corrosion cracking. 11:12 Some 
authors suggested that hydrogen causes embrittlement 
in Pd by decreasing the cohesive strength of GB, since 
H acts electron-receptor and weaken the bonding of host 
atoms across GB. 13,14 

Due to the intricate nature of these phenomena, real¬ 
time measurements and observations can be restricted. 
15 Recent advances in computer powers and computa¬ 
tional theory have made the atomistic simulation an ex¬ 
cellent tool for the investigation of materials properties. 
In fact, first-principles calculations have been widely used 


and succeeded in describing structural properties at elec¬ 
tronic level, particularly for the treatment of defect phe¬ 
nomena in solids. n ’ 16 ~ 21 In this paper, we have exam¬ 
ined the interactions of hydrogen with defects in both 
Pd single crystal and bicrystal. The segregation of H at 
grain boundary and vacancy up to their saturated state 
was investigated in order to understand the process of H 
localization in palladium at the very fundamental level. 


II. COMPUTATIONAL METHOD 

The Density Functional Theory (DFT), 22 as im¬ 
plemented in Vienna Ab-initio Simulation Package 
(VASP), 23 was employed to perform the first-principles 
calculations. In this work, the crystal properties of 
Pd were described by the electron Projector-Augmented 
Wave (PAW) methods 24 with the PBE Generalized Gra¬ 
dient Approximation (GGA) exchange-correlation, 25 
unless noted otherwise. Bulk Pd crystal was simu¬ 
lated using a supercell of 108 face-centered cubic (FCC) 
atoms; in which the convergence of total energy is 0.1 
meV per atom. Additionally, a bicrystal model of 80 
atomic sites was constructed to simulate a symmetric 
tilt £5/(210)/(100) grain boundary in Pd. In the k- 
space sampling, we have used a Monkhorst-Pack grid 26 
of 4x4x4 for the single-crystal and 4x4x2 for bicrys¬ 
tal supercells, for the plane wave basis, respectively. In 
calculating the total energy, equilibrated structures are 
obtained by relaxing the atomic positions and the lattice 
vectors until stress components a xx , a yy and a zz are all 
well below 1.0 kbar. Periodic boundary conditions are 
imposed in three dimensions. 

The precipitation of impurity atoms ensues from the 
migration through grain boundaries, as well as the gen¬ 
eration of vacancies. To demonstrate the binding ten¬ 
dency between hydrogen impurity and palladium host 
atoms, the concept of absorption energy is used; its aver¬ 
age value, F, a b s av9 , from absorption of n hydrogen atoms, 
is calculated by the following equation: 18 
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nKto = E ?o” - El H t -\E H * ( 1 ) 

where E to t nH and E to t 0H describe total energy of systems 
with n hydrogen segregation and at clean state, respec¬ 
tively. E^ 2 is the total energy of one isolated hydrogen 
molecule, calculated by placing that molecule in a large 
supercell (20x20x20 A). 

The formation energy of one vacancy in any system is 
obtained from the change in total energies of that system 
before (E tot ) and after the formation of vacancy (E tot v ), 
adding the average energy of one Pd atom in bulk state 
(E bulk pd ): 

E f = E? ot -E tot + EZ k (2) 

In order to illustrate the site competition between H 
segregation in vacancy (or grain boundary) with H bind¬ 
ing in interstitial, we calculate the trapping energy of 
H. Its average value therefore is calculated with respect 
to the change in total energy of bulk palladium E tot 1H - 
E tot 0H , by adding one H as interstitial: 27 

E v tr ac [nH] = -(E?»- v - E° 0 H t ~ v ) ^ (E} 0 H t - 15°*) (3) 
n 

E tr i n H] = -{E^~ 9b - E? 0 H t ~ sb ) - (E™ - <?) (4) 
n 

where E to t xH ~ v (or E tot xH ~ 9b ) denotes the total energy 
when we have a complex of x H atoms segregated at a 
vacancy (or grain boundary) in the system. The trapping 
energy of the single n th H atom can be estimated by the 
difference in total trapping energies between n and (n-1) 
H atoms: 

A E?r[n] = nE™ c [nH} - (n - 1 )££“[(« - 1 )H\ (5) 
A E*[n] = nE° b [nH] - (n - l)E%{{n - 1 )H\ (6) 

III. RESULTS AND DISCUSSION 

1. Hydrogen binding at interstitial and vacancy 

For a simple crystal such as FCC palladium, octahedral 
(O-site) and tetrahedral (T-site) are two common types 
of interstices. The nearest tetrahedral interstitial (T) is 
in [1 1 1] direction. The nearest octahedral interstitial 
(01) is in [1 0 0] and the second nearest one (02) is in 
[1 1 1] as at the center of the unit cell body (Fig. la). 
Using equation 1, the calculations of bulk Pd show that 
the binding (or absorption) energy of H to O-site is - 
0.1 eV, which is stronger than that at T-site (-0.05 eV). 
This suggests a preference of H occupation at O-site over 
the other in Pd matrix. 28 However, as the magnitude of 
this binding energy is quite low, the stability of H at 
interstitial sites is subject to other factors. 



FIG. 1. (color online) (a) Atomic configuration of a palla¬ 
dium face-centered cubic unit cell, showing tetrahedral (red) 
and octahedral (blue) interstices. In perfect FCC Pd, octa¬ 
hedral site is more energetically favorable for H binding than 
tetrahedral H binding, (b) Variation of trapping energy of H 
along diffusion paths into the vacancy (darker Pd position in 
the figure above). The zero level corresponds to the binding 
energy of H to octahedral interstitial. 


It is believed that hydrogen binding could be stronger 
at other defects such as vacancy and grain boundary. In 
Fig. lb we present the variation of H trapping energy 
along diffusion paths into the vacancy (dark Pd sphere). 
Trapping energy is calculated with respect to binding en¬ 
ergy of one H atom to the octahedral interstitial (the zero 
level corresponds to the binding energy of H to octahe¬ 
dral interstitial). In agreement with previous works in 
metals, 10,27 our calculation also showed that the center 
of vacancy is a very unfavorable binding site for H. Parts 
of the curves below the zero level indicate possible re¬ 
gions around a vacancy at which H binding is stronger 
than that of a general octahedral binding (with no va¬ 
cancy). At the same time, this indicates that the for¬ 
mation of the nearest-neighbor vacancy makes H binding 
at octahedral site stronger (01 in Fig. 1) . However, it 
is interesting that the most energetically stable site for 
H-vacancy binding is in [111] direction, close to its tetra¬ 
hedral position (T in Fig. 1) . This H bonding will be 
threefold, instead of fourfold, due to one missing Pd (va¬ 
cancy). Farther in [111] direction, vacancy has no effect 
on H binding at its second nearest octahedral interstitial 
(as trapping energy at 02 site is zero). 

Using local-density approximation (LDA) method, 
Vekilova et. al 27 reported about a possibility of multiple 
hydrogen occupancy, in which a maximum of six hydro¬ 
gen atoms can be trapped in a monovacancy. The most 
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TABLE I. Formation energy of monovacancy, divacancy and 
influence of hydrogen segregation 


Vacancy type 

Formation energy E/, eV 

Mono-vacancy 

1.17 

Mono-vacancy expr. a 

1.50 

Second vacancy in [100] of VacHo 

1.19 

Second vacancy in [110] of VacHo 

1.13 

Second vacancy in [100] of VacHs 

1.22 

Second vacancy in [110] of VacHs 

1.04 


a Ref. 29 


favorable sites were reported to be along the directions of 
[100] family, as six octahedral sites coordinate a Pd atom. 
However, in our work with GGA, that configuration of 
VacH 6 is not energetically favorable. Instead, we found 
that the stable structure of monovacancy with most H 
will contain eight H atoms, surrounding the vacancy at 
its tetrahedral sites. With this maximum number of H 
occupations, the H binding is still competitively stronger 
than the general octahedral binding in bulk (Fig. 2a). 
However, the insertion of any extra H into this VacHg 
configuration will result in a positive trapping energy, 
i.e. the extra H is likely to diffuse away and occupy an 
interstitial site elsewhere. Fig. 2b shows the variation of 
trapping energy of the 9 th H with respect to its distance 
from vacancy center in [100] direction. There is nowhere 
to fit the 9 th H into this VacH 8 complex to make it a 
stable complex. 

In addition, the vacancy formation process can be 
affected by different factors, including the presence of 
neighboring vacancies and impurity atoms. It is com¬ 
monly believed that the formation of a di-vacancy is 



123456789 
Number of H atoms at vacancy 



FIG. 2. (color online) (a) Variation of H trapping energy as a 
function of H number at vacancy site. A Pd-H configuration 
with hydrogen filling up 8 tetrahedral sites around vacancy 
is energetically favorable and structurally stable, (b) Trap¬ 
ping energy of the 9 th H as function of distance in [100] from 
vacancy center. Its values is everywhere positive along [100] 
direction 


less costly than the formation of two separate vacancies. 
However, our calculations show that this fact is observed 
only in [ 110 ] direction, i.e. di-vacancy formed by the pair¬ 
ing of two nearest neighbors (Table I). The formation of 
the second vacancy in [ 100 ] is shown to bear little to no 
effect. 


2. Hydrogen segregation at grain boundary 

In this subsection, we will discuss the occupation of H 
at a E5 (210) grain boundary, as mentioned earlier that 
grain boundaries play a significant role in the issue of 
impurity transport and segregation. Using equation 4, 
the average trapping energy was calculated as a function 
of H occupation at this GB (Fig. 3). The absorption 
(segregation) energy differs from trapping energy by the 
value of one H octahedral binding. According to Rice- 
Wang theory, 16 , it was commonly accepted that H is an 
interfacial embrittler by segregation energy . 30 

Our simulation showed that the average segregation 
energy is negative for H occupations up to 36 atoms at 
this E5 (210) grain boundary. In other words, a complex 
of GBH 36 could be thermodynamically feasible. How¬ 
ever, Fig. 3 indicated that the average trapping energy is 
positive with approximately more than 28 H occupations, 
, i.e. further absorption of H is energetically unfavorable. 
Fig. 4 shows the atomic configuration of Pd bicrystal cells 
at H-clean state and with 28 H occupation around GB 
plane. If only five Pd atom layers around GB are taken 
into account, then the local H concentration is roughly 
140 % at. (Pd 20 H 28 ) 

Due to the geometry of GB, there are hollow sites on 
GB plane that could host impurity atoms. However, H 
atoms are unlikely to stay in those empty spaces or form 
hydrogen gas themselves; instead, they bind closely with 
Pd. The average H-Pd bond length is approximately 1.8 
A, which is roughly equivalent to a bond between a tetra¬ 
hedral H and Pd in bulk. Also, we found that the most 
common H-Pd bonding at GB is threefold. This con¬ 
figuration contains an H atom bound with Pd at three 



Number of H at grain boundary 

FIG. 3. (color online) Average trapping energy as a function 
of H number at grain boundary region. The absorption energy 
differ from trapping energy by the value of E;, at octahedral 
site in bulk. 
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(a) (b) 



FIG. 4. (color online) Atomic configurations of Pd bicrystals 
at H-clean state (a) and with saturated segregation of 28 H 
atoms. Large spheres demonstrate palladium and small ones 
are hydrogen atoms. Trapping energy of H decreases with 
increase in its local concentration at grain boundary. When H 
local concentration at grain boundary exceeds the saturation 
level, hydrogen is attracted by interstitial sites in bulk. 


vertices of a (quasi-)tetrahedron, while the the forth ver¬ 
tex is missing. This H-Pd binding is identical to the 
configuration of H-vacancy complex, while H found it¬ 
self most comfortable at the tetrahedral interstitial (but 
bonds with only three Pd due to one missing Pd va¬ 
cancy, Fig. la). While H binds threefold, it displays that 
one Pel atom can also participate in several H-Pd bonds 
(Fig. 4b). 

To differentiate between Pd atom layers at grain 
boundary, we index them GBO, GB1, GB2 ..., depend¬ 
ing on how far they are from the grain boundary plane 
(Fig. 4a). Due to the symmetry characterization of this 
GB, layers with opposite indices are equivalent (e.g. GB1 
and GB-1). When H atoms surround Pd, firstly, they 
weaken Pd-Pd bond. Secondly, Pd-Pd can even be bro¬ 
ken, as there is an expansion between Pd layers due to 
high H-GB occupancy. For instance, the distance be¬ 
tween layers GB2 and GB-2 increase from 3.83 A (at El- 
clean state) up to 4.82 A at the suggested H saturation 
(28 H atoms for this 80-Pd bicrystal supercell). 

Next, we investigate the question of how H localization 
would influence the formation of vacancy around the GB 
domain. As presented in Table II, the first Pd layers 
(GB1 and GB-1) are more likely to form vacancy, with 
the formation energy significantly reduced. Farther away, 
the effect of grain boundary and vacancy formation would 
be minimal. The presence of hydrogen in this case, how¬ 
ever, has minimal effect on the formation of vacancy at 
the GB vicinities. For instance, with eight (or twenty- 
eight) H atoms segregation at grain boundary, the va¬ 


cancy formation energy at GB1 even increases slightly to 


TABLE II. Formation energy of monovacancy at different 
sites around grain boundary, eV 


Vacancy site 

Formation energy E/, eV 

GBO 

1.29 

GB1 

0.48 

GB2 

1.00 

GB3 

1.04 

GB4 

1.06 

GB5 

1.09 


0.55 eV (or 0.63 eV, repectively), compared to its value 
of 0.48 eV in H-clean state. Therefore, unlike H-vacancy 
binding, H-GB segregation doesn’t tend to initiate hy¬ 
drogen embrittlement through microvoid formation. In¬ 
stead, high H content could result in a decohesion across 
the Pd grain boundary and subsequently, the hydrogen 
embrittlement could be initiated by the grain separation 
at sufficient external stress. This observation is in agree¬ 
ment with calculations using Rice-Wang model, in which 
single H was reportedly a GB embrittler by reducing the 
grain cohesion. 30-32 

IV. CONCLUSIONS 

We have reported the saturated localization of H bind¬ 
ing at vacancy and grain boundary in Pd crystals. The¬ 
oretical investigations for favorable configurations of H- 
vacancy and H-GB complexes from low to high H occu¬ 
pation were presented. 

While H prefers to occupy the octahedral interstitial 
position, they are anticipated to fill in the tetrahedral 
sites in complex with a Pd vacancy. Our calculation 
showed that one Pd vacancy is capable of storing up to 
eight H, instead of six as reported previously. 27 Pairing of 
vacancies is encouraging in [110] direction, and H pres¬ 
ence will further assist multi-vacancy formation, which 
implies the hydrogen embrittlement by connections of mi- 
crovoids. Numerous experimental observation and theo¬ 
retical calculations have reached to the same conclusions 
on vacancy formation induced by hydrogen. 

In both cases of open defects (vacancy and grain 
boundary), hydrogen prefers a threefold bonding with 
palladium. However, H segregation at GB results in 
weakening Pd-Pd bonds across GB, by inducing the grain 
separation. This observation supports the decohesion 
mechanism of hydrogen embrittlement as evidenced by 
many experiments. 
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